function     PINAW = calc_PINAW(actual, lower, upper)
    interval_widths = upper - lower;
    PINAW = sum(interval_widths) / (length(actual) * (max(actual) - min(actual)));
    % PINAW = mean(upper - lower);
end